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Abstract: By studying the Euclidean partition function on a cone, we argue that pure 
and mixed gravitational anomalies generate a "Casimir momentum" which manifests itself 
as parity violating coefficients in the hydrodynamic stress tensor and charge current. The 
coefficients generated by these anomalies enter at a lower order in the hydrodynamic gradient 
expansion than would be naively expected. In 1 + 1 dimensions, the gravitational anomaly 
affects coefficients at zeroth order in the gradient expansion. The mixed anomaly in 3 + 1 
dimensions controls the value of coefficients at first order in the gradient expansion. 
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1. Introduction 

Anomalies constitute a set of fascinating phenomena in field theory. Their usefulness stems 
mainly from their robustness — they are often easily computed even in very strongly interacting 
theories and they are insensitive to a large variety of deformations of the theory. Further, 
their phenomenological significance in fields ranging from solid state physics to cosmology has 
inspired a large body of work in the last few decades devoted to understanding the dynamics 
of anomalies from various viewpoints. 

Despite this effort there is much that is yet to be understood. We have but a vague 
comprehension of the dynamics of anomalies at finite temperature and chemical potential, 
let alone far away from equilibrium. The existence of thermodynamic observables which are 
sensitive only to anomalies, or a classification of hydrodynamic processes which are responsive 
only to anomalies are examples of problems which one would like to address. 

While a complete understanding of the role of anomalies in terms of thermodynamic and 
hydrodynamic response is still elusive, a significant amount of progress has been made in the 
last few years using only a handful of tools and techniques. These tools have successfully 
explained various features of anomaly induced response in the case of pure U(l) anomalies 
in hydrodynamics [1]. In what follows, we would like to argue that the existing methods 
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which have been so successful in dealing with U(l) anomalies are bound to give incomplete 
information when applied to gravitational or mixed anomalies. In order to do so, let us briefly 
summarize some of the salient features of hydrodynamics. 

The equations of motion of relativistic hydrodynamics are energy-momentum conserva- 
tion and charge conservation when appropriate. These equations must be supplemented by 
constitutive relations which express the (covariant) energy momentum tensor T^ u and current 

in terms of the hydrodynamic variables: T, fi and u M — temperature, chemical potential 
and velocity field respectively. We normalize the velocity field such that u 2 = — 1. If we also 
turn on a slowly varying background gauge field and metric g^, then the constitutive 
relations will depend on gauge and diffeomorphism invariant combinations of the background 
fields. 

The constitutive relations can be decomposed into scalar, vector, and tensor relations 
with respect to the rotational symmetry which preserves the velocity field, 

= Nu» + v» , 

(1.1) 

= £u»u v + Vb!™ + u»q v + u V 

where we have defined the transverse projector = g^ u + u^u v . The vectors and 
are transverse to u^, and t^ v is transverse and traceless. The constitutive relations (|1.1|) are 
usually arranged in a derivative expansion in which expressions which involve gradients of the 
hydrodynamic variables are suppressed by appropriate powers of the mean free path [2,3]. 
The constitutive relations characterize the response of the energy momentum tensor (or of 
the charge current) to a perturbation via response parameters. For instance, the conductivity 
parameterizes the response of the charge current to an external electric field. 

In some instances "transport coefficients" and "response parameters" are used inter- 
changeably. In what follows we will make a distinction between transport coefficients or 
hydrodynamic response parameters, and thermodynamic response parameters (whereby re- 
sponse parameters or response coefficients will refer to either set). Transport coefficients 
characterize out of equilibrium dynamics and their Kubo formulae will involve correlators at 
non zero frequency. Thermodynamic response parameters are associated with the properties 
of the equilibrated system and the the appropriate Kubo formulae will involve correlators at 
zero frequency. The shear viscosity, bulk viscosity, and conductivity are examples of transport 
coefficients. An example of a thermodynamic response parameter is the magnetic suscepti- 
bility. 

We are now in a position to discuss anomaly-induced transport as it is currently under- 
stood [1,4-22]. In 1 + 1 dimensions, we parameterize the anomalies through the conservation 
equations for the covariant current and stress tensor, 

where F^ u is the field strength of the background gauge field, R is the Ricci scalar of the 
background metric, and V indicates the covariant derivative using the Levi-Civita connection. 
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To zeroth order in derivatives, the constitutive relations are given by 
V = P, £ 

■ at 

v» = -2c s ^ u u v , tt = e^u u (c 2d T 2 



dP BP 



C s fJ, 



0. 



(1.3a) 
(1.3b) 



where P is the pressure and t^ v identically vanishes in two dimensions [15,21]. The param- 
eter C2d is a free constant which violates parity and time-reversal while preserving charge 
conjugation. At the one and two derivative level equation (|1.3|) receives several corrections 



which we will discuss in section 4.1. The author of [23] showed that in the presence of a two 
dimensional gravitational anomaly the vector contribution q 11 to the constitutive relations 
( |L3bD b 

/ \— X~t>T\ 

(1.4) 



In 3 + 1 dimensions we parameterize a U(l) 3 and mixed gauge- gravitational anomaly by 

1 



^CAF^Fpcr + c m R a /3lll/ R apa 



(1.5) 



To first order in derivatives, the constitutive relations are given by [1, 10, 11] 



v = p- cvx , 



dP BP 



dP 



(1.6a) 



and 



^ = -QcaiiB^ + (c 4d T 2 - Scau 2 )^ + <jA^ ( E v - TV V 
q» = (c4dT 2 - 3c^ 2 )5" + 2(c 4df iT 2 - c^K , 



(1.6b) 



-7]<J 



where we have defined 

£M = ^P- Uu F pa , 



E„ 



F ii u 



a 



e» v P°u v V p u a , 



(1.6c) 
(1.6d) 



From the hydrodynamic analysis of [10], the parameter c 4d is a constant which violates parity 
and charge conjugation while respecting time-reversal. There is no four dimensional analogue 



x The gravitational anomaly coefficient c g is denoted by the symbol —D in [23]. Further we have used the 
following equilibrium relation to simplify the result derived in [23] 



V 2 T 
T 
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of the two dimensional computation of [23] in the literature. We carry out such an analysis 
in this paper and our results are presented in appendix ||. 

Explicit calculations in various 1 + 1 dimensional CFTs with pure gravitational anomalies 
support the relation 

C2d = -8TT 2 Cg . (1.7) 

The evidence for this statement in various 1 + 1 dimensional CFTs is extensive, though the 
results are somewhat fragmented across papers in various subfields of physics, ranging from 
thermal hall physics and superfluids (See for example [24-27] and references therein) through 
free field computations [28,29] to various studies of AdSs/CFT2 and thermodynamics of BTZ 
blackholes in the presence of gravitational Chern-Simons terms [30] (See [31] for a review). 
Similarly, computations in an arbitrary free field theory of chiral fermions [26, 29, 32-37] (see 
also [24]) and in AdS5/CFT4 [38] are consistent with the relation 

C4d = -87r 2 c m , (1.8) 



in four dimensions. A natural generalisation of ( |l.7D and (1.8) to arbitrary dimensions was 
presented in [29] by an analysis of free fermion theories. 



The peculiar relations (|bj) and (|1.8j ) are the main focus of this paper. They exhibit two 
distinct features which a priori indicate that a general proof of their validity (if it exists) will 
require special attention. 

1. A breakdown of the naive derivative expansion: The gravitational and mixed anomalies 
contribute to the response parameters at two orders of derivatives lower than expected 
when studying the equations of motion. We will refer to this phenomenon as a "jump" 
in the derivative expansion. For instance, the gravitational anomaly in 1 + 1 dimensions 
enters at zeroth order in the derivative expansion instead of second order in the deriva- 
tive expansion. Existing methods [1, 19, 39] independently constrain the constitutive 
relations order by order in derivatives, and so cannot exhibit jumps in the derivative 
expansion. Indeed, when applied to the case of a two dimensional theory with a gravita- 
tional anomaly [23], existing methods lead to the contribution ( |1.4[ ) to the constitutive 



relations, but fail to establish (1.7). 



2. The 8ir 2 factor: All existing methods to derive constraints on anomaly induced trans- 
port [1,10,11,19,39] carry out algebraic manipulations on the anomalous conservation 
equations (for example ( |1.2|) ). Using our normalization for the various coefficients, a 
moment of thought reveals that the constraints derived out of such algebraic arguments 
cannot contain transcendental numbers like the additional factor of 8tt 2 in ( |1.7j ). This 
obstruction by itself precludes any constraint of the form (|0]) to arise from methods 
based on [1]. Moreover, any argument which gives a constraint such as ( |1.7[ ) or ( |1.8| ) is 
expected to be geometric rather than algebraic. 

Fortunately, there is a very well-known relation in 1 + 1 dimensional conformal field theory 
(2d CFT) which exhibits the same features as ( |1.7|) or ( Oj ) — the Cardy formula. Given a 
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unitary modular-invariant 2d CFT with a Weyl anomaly 

(T cmj )» = c w R=±^^R, (1.9) 

the Cardy formula relates the pressure of this CFT in the infinite volume limit to the left and 
right handed central charges cr and cl via [40,41] 

P = 2nT 2 ^^ = Air 2 c w T 2 . (1.10) 

The relation ( 1.1CD follows, essentially, from the Casimir energy picked up when going from 



two-dimensional Euclidean space to the cylinder due to the non-conformal behavior of the 
energy momentum tensor. The relation ( 1.10| ) exemplifies a relation between terms of different 



orders in the derivative expansion. (In this case, the pressure is zeroth order in the derivative 
expansion while the Weyl anomaly is second order.) Further, the transcendental number here 
47T 2 comes essentially out of the geometry of the thermal circle; 2-kT is precisely the inverse 
radius of the cylinder. In fact, as we show in appendix [A], one can easily extend the result 
Ql-lCp to theories with c g 7^ therey deriving (1.7) for conformal field theories with no charge 



current. 

In the rest of this work, we use the intuition gained by the observation ( |1.10 ) to demon- 



strate (P-.7D and ( |1,£| ) by generalizing some recently obtained results for the Euclidean gener- 
ating functional [19,39]. We will argue that the jumps in the derivative expansion in 1 + 1 
and 3 + 1 dimensional theories can be understood in terms of a "Casimir momentum density" 
which can be computed by placing the theory on a cone (or a product of a cone with a two 
dimensional manifold). We also compute the constitutive relations of a conformal and non- 
conformal theory in 1 + 1 dimensions to second order in the derivative expansion (see ( |4.7| )), as 
well as the contribution of the mixed anomaly to the constitutive relations in 3 + 1 dimensions 
(which may be found in appendix [B]). 

The remainder of this manuscript is organized as follows. In sections || and || we collect 
various results from the literature relevant to this work: section |2| briefly reviews the structure 
of the anomaly polynomial and the construction of the Bardeen-Zumino anomaly polynomial 
which distinguishes between covariant and consistent observables. Section || summarizes 
the method of construction of an equilibrium partition function [19, 39] which will play an 
important role in our work. The reader familiar with the content of sections ||-|3] may go 
directly to sections ^| and || which contain our main result. In section |] we compute the 
effect of a gravitational anomaly on the response coefficients of two dimensional theories and 
in section |5| we compute the effect of a mixed anomaly on the response coefficients of four 
dimensional theories. We end with a discussion and summary in section ||. 

Note: While this work was nearing completion, we were made aware of [42] which deals 
with issues similar to the ones considered in this paper. 

2. Anomalies 

In this section we briefly review several useful facts about gauge and gravitational anomalies 
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and establish our notation. Since these facts are well known, we will necessarily be brief. We 
will refer the reader to various books and lecture notes [43-45] for a more detailed review. 
The reader familiar with these results may skip to Section |3[ 

2.1 Anomaly polynomial, Consistent and covariant anomalies 

Consider a 2n-dimensional theory with a conserved U(l) current and stress tensor which 
respectively couple to a background gauge field A and metric g. We will assume that these 
conservation laws become anomalous once we place the theory on a non-trivial background. 

We will denote the curvature of A by F = dA, and we define the Christoffel connection 
1-forms and curvature 2-forms associated with the metric via 

r^rv^, K a (s = \R a ^ v dx» a dx v = dv a p + r% a r>, (2.1) 

where Greek indices are spacetime indices, are the Levi-Civita connection coefficients, 

and R a pav are the components of the Riemann tensor. 

Let W denote the generating function of the 2n dimensional theory we are interested in. 
Let us also define a 2n + 1 dimensional theory on a manifold A^n+l with boundary dM.2n+i 
with a generating function W cov defined through 

W cov = W[dM 2n+ i] + / I 2 C n+i • (2-2) 

J M'2n + 1 

Here I^n+i ls a Chern-Simons form associated with the metric, gauge field, or a combination of 
the two. The form of the 2n + 1 dimensional Chern-Simons term is determined by demanding 
that the covariant generating function, W cov , is invariant under gauge transformations and 
diffeomorphisms. More precisely if we denote gauge and diffeomorphism variation by 5\, 

5 X A^ = d^A + A u d„e + {d v A„)e = d, L (A + A v ?) + 

bxg^u = + , 

then the statement that the 2n + 1 dimensional theory is gauge and diffeomorphism invariant 
implies that 5\W cov = 0. 

As a result of the gauge invariance of W cov , the 2n dimensional theory described by 
the generating function W in ( |2.2[ ) is not gauge and diffeomorphism invariant and hence, 
anomalous. In detail, the gauge and diffeomorphism invariance of W cov amounts to 

5 x W[dM 2n+ i\ + h [ lg S +1 = . (2.4) 

By construction, the bulk Chern-Simons term is gauge-invariant up to boundary contribu- 
tions. Let us parameterize the variation of the Chern-Simons term by 

Sx [ lg S +1 = [ d 2n x^~g A J + / d 2n x^~g d x C T\ , (2.5) 
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so that (2.4) reads 

5xW[M 2n ] 



J 

J Mi 



d 2n x^~g A J 



J 

J M2 



d 2n xj=g~ dxC T\ + 0. 



(2.6) 



The fact that the generating function W is not gauge and (or) diffeomorphism invariant 
signals the presence of anomalies. From the viewpoint of the 2n + 1 dimensional theory, the 
2n dimensional anomaly is due to an anomaly inflow — a flow of conserved charges from the 
bulk to the boundary thereby modifying the conservation laws at the boundary. 

Let us consider these conservation laws in some detail. We define a consistent stress 
tensor and a consistent charge current by varying the generating function W with respect to 
the metric and gauge field respectively. 2 



u cons 



rpflV 



1 SW 

2 5W 



(2.7) 



We can now compute the gauge and diffeomorphism variation of W, 



SxW 



d 2n x 



d 2n x. 



SW . SW . ' 
OxQfiu + jj-oxAy, 



{21 



g [ A(V^ Jcons) £tf K^v^oons ^^vJ C ons + A^V uJ C ons)] 



where in going from the first line to the second we have integrated by parts. Here V represents 
the covariant differentiation of tensors using the usual Christoffel connection. Generically, in 
a theory with diffeomorphism anomalies none of the currents are tensors — we define the action 
of V on Tco ns and Jcons by treating them as if they were tensorial objects. Comparing ( |2.8| ) 
against (P]q) we get the consistent anomaly equations: 



V u Jlt 



n ° cons 



J, 



V7 rpnv _ -pn TV 



v" cons 



JA^-g 



,flU 



-9 



dx 



-gT> 



(2. 



We emphasize that neither the currents Jcons and Tcons nor the anomalies on the right 
hand side of equation ( |2.9| ) transform covariantly under gauge transformations and diffeomor- 
phisms. This unfortunate state of affairs may be remedied by using the generating function 
Wcov instead of W. 

Consider a general variation of the 2n + 1 dimensional metric g^ +1 ^ and the 2n + 1 
dimensional gauge field A^ n+1 ^ on M.2n+i- I n what follows Roman indices will be associated 
with the coordinates on M.2n+i and Greek indices with the coordinates on dM.2n+i- The 



The "consistency" of the current and stress tensor are associated with the so called Wess-Zumino consis- 
tency conditions [46] (See [43-45] for a review). 
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index associated with the direction normal to the boundary will be denoted by JL. The 
orientation of the normal direction is fixed by demanding 



- 5 (2n+l) e J-W- i -M2n = ^ e Mi-M2n _ ( 2 .10) 

The variation of the 2n + 1 dimensional metric will contain a bulk piece bg ab and a boundary 
piece bg^ v . Likewise, the variation of the gauge field will contain a bulk piece SA a and a 
boundary piece SA^. With this notation in mind, we can parameterize the variation of W cov 
with respect to the sources through 

b S [ I 2 C „+i = I d 2n+l x^-g^)h a 5A a + L ab c ST c ab \ 

J M 2n + 1 JM 2n +l " ^ J 

+ / d ^ x ^\ Y ^A^ + lx x ^5r\A. (2.11) 

JdM 2n + l 1 1 } 

Before we proceed, certain clarifying remarks are in order. The metric variation enters 
entirely through the variation of the connection because of the topological nature of the Chern- 
Simons terms. Furthermore, we have ignored all dependencies on, say, the extrinsic curvature 
of the boundary since it is not relevant to the discussion below. 3 Since the Christoffel symbols 
are symmetric in their last two indices, the tensors L and X are well defined only up to tensors 
which are anti-symmetric in their first two indices. In what follows we will fix this ambiguity 
by demanding that 

Labc = -L ac f>, . . 



Converting the variation of the Christoffel symbols in (2.11) to metric variations we can 
write 

bs I I& s +1 = I d 2n+1 x^g^-Ar5A a + V c L abc 6g ab \ 

JM 2 n + l JM 2n + l 1 J 

+ / d^x^-gip^bA^^-P^bg^X (2.13) 

JdM 2n+ i 1 1 ) 



with 



P» Z = Y*, P^^-^Vx^X^ + X^-X^, (2.14) 

and Greek indices are raised and lowered with the metric g^ v . In the equations above we 
have again dropped extrinsic contributions to P^ z (of the form + L u±fl coming from an 
integration by parts) . 



3 This involves the following assumption. The contractions in the boundary terms involve only the sum over 
the boundary indices. While this is sufficient for our purposes, we warn the reader that a more careful treatment 
of such extrinsic terms is required in many other instances. For example when carrying out computations in 
the framework of the AdS/CFT correspondence such terms (along with other Gibbons-Hawking like terms) 
give a finite gauge(or diffeomorphism)-invariant contribution to the dual stress tensor [30]. 
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We are now in a position to define the covariant currents as the variation of W a 
respect to the gauge fields and metric on dM.2n+ii 



with 



_ 1 5W cov 



J- COV 



2 5W CI 
^ Sg. 



TM _i_ pt 1 
"cons ' 1 BZ' 



A cons ' r BZ ' 



(2.15) 



[IV 



By construction these currents transform covariantly under gauge transformations and dif- 
feormorphisms (they were obtained by varying an invariant functional). In other words, the 
non-covariant transformations of Pjg Z and P i B w z under gauge transformations and diffeomor- 
phisms exactly compensates for the non covariant transformation properties of the consistent 
currents, implying that the covariant currents transform covariantly as advertised. The cor- 
rections P BZ an d P&z that come from varying the bulk Chern-Simons terms are are often 
termed the Bardeen-Zumino currents or Bardeen-Zumino polynomials [47] which is the reason 
for the subscript l BZ' in our notation. 

Let us turn to the conservation laws for the covariant currents. To compute these con- 
servation laws we need to express the divergence of the Bardeen-Zumino currents in terms of 
T^ u and J defined in (2.5). This computation can be carried out by using ( 2.13| ) to evaluate 
the gauge and diffeomorphism variation of f I^n+i- Since the variation of the Chern-Simons 
term gets contributions only from the boundary we find, 



M 2n +i 



T cs 

'271+1 



dM 



2n+l 



/ ^W^AtV^-J^. (2.16) 



Comparing ( |2.16[) against (^) we get 



V pi 1 — t-L _ j 



V v Pg z = F%P U BZ + JA» + V U L L ^ + g 



0^ 



~gT> 



(2.17) 



Finally, combining ( |2.17 ) with ( |2.9D we find 



V7 TV — t-L 

v fJ. u cov ~ ° 1 

V7 j<(iv _ z?M f _|_ v J^tiv 

V V- L COV ~ 1 V^COV ' V v i-l 



(2.18) 



An alternate way of deriving the conservation law ( p. 18 ) is to begin with the variation of 
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Wcov with respect to an arbitrary variation in sources 

6sW cov = 6 s W[dM 2 n+i] + 5s / ^2n+l 

I d 2n+1 x v ^^7{™ a + V c L afec ,5 5af) ) 

/ d 2n x^-g{j£ ov 5A a + l -TZ5guv\i (2.19) 

JdM 2 n+i 1 1 J 



IM2n + l 

+ 



and then demand that the gauge and diffeomorphism variation of W cov vanish. From this 
derivation it is clear that the covariant anomaly equations ( |2.18| ) can be thought of as the 
conservation laws for the covariant charge applied to M.2n+i and its boundary. This anomaly 
inflow phenomenon is often called the Callan-Harvey anomaly inflow mechanism [44,48-50]. 

While not relevant for this work, before we move on to particular examples, we note that 
we can always define a canonically conserved energy-momentum tensor 

rpfM _ rpUV T -LflV 

conserved — 1 cov ^ J 

such that 

V JC — T-L 

v ii° cov ° ) 

^ 7 ^conserved = ^^f^cov^ (2.20) 

rpflll rpUfJ, nT^flV 

conserved conserved ' 

where the last equation can be interpreted as the Lorentz anomaly — the non-conservation 
of the angular momentum current. Thus, —2L- LfJ-u can be thought of as the rate at which 
angular momentum is injected into the system. This demonstrates the well-known result in 
the theory of gravitational anomalies: we can always trade a diffeomorphism anomaly for a 
Lorentz anomaly. 

Note that to study gravitational anomalies one needs to necessarily place the 2n dimen- 
sional theory on a non-Minkowski background without Poincare invariance and hence we need 
to clarify what we mean by 'angular momentum' in the discussion above. Consider a Killing 
vector which generates a symmetry of the given gauge and gravitational background. More 
precisely, we define a vector and a corresponding U(l) gauge transformation parameter 
Ak such that 

£k9uu = ^uK v + V U K U = 0, 
£ K A li + V IM A K = K v F vli + V ll (K v A 1 , + A K ) = Q. <2 '" ' 

According to Noethers theorem there exists a Noether Current 

J£ = K u TZ served + {K V A U + A K ) J? ov , (2.22) 
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whose diveregence can be worked out using ( 2.20 ) to be 

1 

2 



V M J £ = o T cLserved £ K9^u + (£kA^ + V^A^) «/£,, 



(2.23) 

-L^(V M ^ - VvK^) + J x (K u A v + Ax). 



The first line of the above equation vanishes because of our assumption in ( 2.21 ) whereas 
for the second line to vanish either the background should have the anomalies turned off 
i.e., 3 ± = and Xj^ v = or an extra restriction needs to be satisfied by A^ and K^, i.e., 
V„K V - V V K„ = and K V A V + A K = 0. 

Thus, the Noether currents corresponding to curl-free Killing vectors ('momentum-like 
Noether charges') can still be conserved in the presence of gravitational anomalies whereas the 
Noether currents corresponding to Killing vectors which are not curl- free ('angular momentum- 
like Noether charges') get a contribution to their divergence proportional to L J -^ tl/ . This is 
the precise meaning of the statement that the theory placed on a non-trivial gauge and grav- 
itational background exhibits a Lorentz anomaly given by — 2L" L/W . 

So far, we have characterized gauge and diffeomorphism anomalies in 2n dimensions using 
a 2n + 1 dimensional Chern Simons form. The most convenient way to characterize 2n + 1 
dimensional Chern-Simons forms is to construct a 2n + 2 form 

Vanom = dlg S +1 . (2.24) 

which we call the anomaly polynomial of the theory. This anomaly polynomial is an index-like 
object built out of the Chern classes of F and Pontryagin classes of 91. Given an anomaly 
polynomial V a nom encoding the anomalies of a theory, one can reconstruct the Chern-Simons 
form an d from it the anomalous behavior of the stress tensor and current. 

2.2 Two dimensional theories 

The anomaly polynomial of an arbitrary theory in d = 2 is a 4-form given by 

Vanom[F, £H] = c s F A F - 8TT 2 c gPl (JR) = c s F AF + c 5 tr(9* A SH), (2.25) 
where p\ is the first Pontryagin class 

Pl (K) = --^tr(WA!K) J (2.26) 

and tr denotes a trace over the Lorentz indices. The coefficient of the pure U(l) 2 Schwinger 
anomaly is denoted c s and c g is the coefficient associated with the pure gravitational anomaly. 

Given an arbitrary theory with chiral bosons and (or) fermions, c s and c g in the anomaly 
polynomial can be computed using 

° s = ~2'(2vr) 2 ^ Xi(qf) 1/2 , 

^ l=SpeaeS (2.27) 

i=species 
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where \i denotes the chirality (we assign right handed fermions positive chirality), q denotes 
the charge and the subscripts {0, 1 /2} represent the contribution from a chiral scalar and a 
Weyl-fermion respectively. The contribution of a Major ana- Weyl fermion is half that of a 
Weyl-fermion of the same chirality. The sum above is performed over each species where each 
particle/anti-particle pair contributes one term to the sum. In particular for a 2d-CFT with 
a £7(1) symmetry, we have 



2vr 
'4!(8^2 



(k R - k L ), 
(cr -c L ), 



(2.28) 



where k R / L are the right/left £7(1) Kac-Moody levels and cr/l are the right/left Virasoro 
central charges. 

The Chern-Simons form corresponding to the anomaly polynomial (2.25) is 



I% s = c s Af\F + c g jcs- 
where we have defined the gravitational Chern-Simons 3-form 



(2.29) 



3CS = tr 



rAdr + -rATAr 

3 



(2.30) 



such that djcs = tr(9t A !SR). The gauge and diffeomorphism transformation of these Chern- 
Simons terms is given by 



M 3 



T CS 
J 3 



dM 3 L 



c 3 AF + c g (d x e)dT X u 



Comparing (2.31) with (|2.5|) we get 



(2.31) 



J = c s ~e a ?F a p , 
T\ = c a e^d n T } 



(2.32) 



Thus, 



V — r -f a @ F a 



Vrp^iv TPf MU jcons 



c s -e^F^AV - Cg gr-j= d. 
2 v 9 



A 



-9* 



u a L vt 



(2.33) 



To obtain the covariant current we vary the Chern-Simons forms with respect to the 
sources, 

5 S [ I$ s = f \2c s 5AAF + 2c g tr(5T AfR)\ + f \c s 5A A A + c g 5T" x A r A „) . (2.34) 
J Ms J Ms *- J JdMs J 
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Comparing fl2.34|) against ( 2.11 ) we can read off the covariant anomaly and the Bardeen- 
Zumino currents 



1 



J- = 2c s -e at5 F al3 , 
h^ v = 2c g ^R^ a0 , 

r BZ — s u > 

Y^ X — o f Wr J -1- r- f A T' 1 
yv v — w/ c y up ~ L j c L upi 



p pu 



_y ( j^Xpu _|_ j^Xup j^puX 



(2.35) 



The covariant anomaly equations (2.15) become 

V 7^ — 2r -f a/3 F a 

v P J cov ~ 2 r aj3i 



(2.36) 



In two dimensions the Riemann tensor satisfies an identity 



which can be used to bring the above equations to a more familiar form 

V M J^ = 2c s e a(3 d a Ap, 

y rppu _ rpp jp , uuy n 



(2.37) 



2.3 Four dimensional theories 

We now turn to four-dimensional theories. The anomaly polynomial 6-form of an arbitrary 
4-dimensional theory can be written in the form 

'Panom[F } = c A F A F A F — 8ir 2 c m F A p x = c A F A F A F + c m F A tr(9l A 91), (2.38) 

where c A is the U(l) 3 triangle anomaly coefficient and c m is the mixed [/(l)-gravitational 
anomaly coefficient. Given a theory with chiral fermions, these coefficients can be calculated 
via 



2tt 



3!(2vr) 3 . 



Yl Kit 



f)l/2) 



i=species 



2tt 



4!(87r 2 )(27r) . 



Xi(«)l/2 • 



(2.39) 



i=speaes 

The Chern-Simons form corresponding to the anomaly polynomial above is 



T cs 

2 5 



A A [c A F A F + (1 - a)c m tr(V\ A JR)] + ac m F A j CS - 



(2.40) 
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where jcs is the gravitional Chern-Simons 3-form defined in ( |2.30D and a is a parameter 
which determines how the mixed anomaly is shared between the £7(1) and the gravitational 
transformations. It corresponds to a gauge and diffeomorphism non-invariant contact term 
in the consistent generating function, W a = — J ac m A A jcs- At a = 1 the mixed anomaly 
is completely associated with diffeomorphism transformations whereas for a = the mixed 
anomaly is completely associated with £7(1) gauge transformations. This becomes clear if one 
works out the gauge and diffeomorphism transformation of the Chern-Simons form 



f lg s = [ A[c A FAF + (l-a)c m tr(!RASR)]+ / a c^F A dT\(d x C) - (2.41) 
JM 5 JdM 5 JdMs 



Comparing (2.41) with fl2.5|) we get 



ij _ Koa.fi 

J ~ 4 



c A F K (jF a p + (1 — a)c m R u \ KU R ua p 



(2.42) 



so that the consistent anomaly equations (|2.9|) become 



v A 1 cons j c 

V7 t^ v = l v 



c A F K aF a fi + (1 — a)c m R u \ Ka R ua p 



c A F Kr jF a /3 + (1 — a)cmR v Xk<jR va.fi 



A* 



a c m g 



au 



Ox 



\e K ^F Kn d n T 



Ko-fJaL vfi 



(2.43) 



The variation of the Chern-Simons term with respect to sources is given by 



6 S 



[ lf 5 = / \6AA[3c A F AF + c m tr(<RA*R)]+2c m tr{5T AOS) AF\ 



+ / 5A A {2c A A AF + a c m jcs] 

JdMs 1 J 



+ 



dM; 



5T U X A \a c m F A T x v + 2(1 - a)c m A A 



(2.44) 



from which we can read off the covariant anomaly and the Bardeen-Zumino currents using 
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4 £ 



3c A F Ka -F a i3 + c m R \ Ka R poifi 



1 



Jos 



pfiu 



c fJ,UK(T 



a c r , 



-iA "pet 



t 1 I/O T c 1 



"P 

1 



+ 2(1 - a)c m - [e^^^ + e A ^ CT i^, 



1. 



V A ( X A ^ + - X^ A ] . 



The covariant anomaly equations ( |2.18| ) are now given by 

1 _ 



V 7^ — -f 

v H u cov — ^ e 



F Ka F a j3 -\- C m R Xko-R ua(3 
1 



a/3 



(2.45) 



(2.46) 



Note that the final covariant anomaly equations do not depend on the parameter a, which 
should not surprise us as a corresponds to a local counterterm. We see that the mixed 
anomaly is symmetrically shared between the U (1) and gravitational currents. 



3. The Euclidean generating functional 

In this section we review the method by which the Euclidean partition function for theories 
in thermodynamic equilibrium may be constructed in terms of a functional of the sources. 
Most of the material contained here is a shortened version of [19, 39] to which we refer the 
reader for more details. 

3.1 The generating function 

We begin by reviewing the arguments presented in [39] used to construct the generating 
function of a thermodynamic theory. Consider a Lorentz-invariant quantum field theory in 
d space-time dimensions at nonzero temperature T in flat space. Generically, the real-space 
Euclidean correlators of this theory will decay exponentially at large distances implying that 
the screening lengths & of the theory, 



(O(r,x)O(r,0)>~exp(-|x|/0, 



(3.1) 



are finite. The assumption that the correlation length is finite implies that, for example, we 
are not at a critical point of the theory, or that there are no unscreened long-range forces as in 
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QED. Instead of screening lengths, one may consider screening masses which are given by the 
location of the poles of the momentum-space zero-frequency Euclidean correlation functions 
along the imaginary momentum axis. The screening lengths are inversely proportional to the 
screening masses. 



From (3A.) it follows that the zero- frequency correlation functions of the theory are an- 
alytic at zero spatial momentum. We then define truncated correlation functions by Taylor- 
expanding all zero-frequency correlation functions to m th order around zero spatial momen- 
tum. The position-space truncated correlators may then be obtained by varying a generating 
function: an integral of a local functional of background fields, which we denote W m . In a 
theory with local gauge and diffeomorphism invariance, W m may be written as an integral 
of gauge and differmorphism invariant (up to boundary terms) scalars. In the presence of 
anomalies, appropriate local (but Lorentz breaking) terms may be introduced into W m in 



order to account for its anomalous variation as described in section 2.1 



Due to the non-analytic behavior of Euclidean correlators in the complex momentum 
plane, one can expect the derivative expansion to have a finite radius of convergence at 
best. Nevertheless, one may always carry out a formal expansion of the Euclidean generating 
function around the origin to arbitrary order. A resummation of this series will presumably 
uncover the poles or branch cuts of W. Terms which are non-analytic at zero momentum, 
e.g., e -1 / fc , will not be accounted for in our construction and must be included by hand when 
studying properties of the theory which go beyond the derivative expansion. We will discuss 
the role that such terms play in our analysis in sections |] and ||. 

Put differently, were we able to resum the derivative expansion and uncover all non- 
analytic contributions, we would have obtained the equilbirum partition function. By equilib- 
rium partition function we mean the Euclidean partition function evaluated on a background 
with time independent sources. In practice this resummation is unfeasible, but even in con- 
figurations where the derivative expansion breaks down one may still extract information 
about correlation functions without performing the resummation. We follow this method in 
sections | and |. 

To make the construction of W m explicit, consider a Euclidean theory defined on a 
manifold Md = S 1 Xf -Md-i, where the S 1 is the time circle which may be fibered non- 
trivially over the spatial manifold Ai^-i- Let t denote time and V = dt, i.e., is the 
vector whose integral curves give the time circle. In thermal equilibrium we demand that 
£y = 0: the Lie derivative of all sources with respect to V vanish. Consequentially V is a 
Killing vector. The generating function W m may be constructed of all possible gauge and 
diffeomorphism invariant combinations of the metric, gauge field, Killing vector V and non- 
local quantities which may be constructed using V, namely the inverse length of the time 
circle T~ l = \/—V 2 dr, which we identify with the temperature, and the Polyakov loop 
Pa = exp(J AfjY^dr). Here r = it denotes Euclidean time. We find it convenient to replace 

with the fluid velocity = V^/yZ—V 2 and the Polyakov loop with a chemical potential 
H = Tin Pa. 
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Our conventions are such that varying W with respect to the sources as in fl2.7| ) will result 
in real one-point functions in Lorentzian signature. That is, we work with a Wick rotated 
generating function. If we denote the partition function of the theory by Z then we define 
W = — i In Z. By construction, all the covariant n point functions transform covariantly when 
Wick rotating to or from Euclidean time. 

Let us work out a simple example in detail. We wish to construct the generating function 
of a parity preserving theory to first order in the derivative expansion. To do so, we must 
enumerate all possible scalar expressions which have zero or one derivatives acting on gauge 
and diffeomorphism invariant combinations of T, /x, g pu and A p . Consider the identities, 

V^iiy = -u^ + tt^u, V (U T = -Ta M , V^fi = -fia fl + E^, (3.2) 

where we have defined 

a M = «"V I/ « M , = ^A^A^(V p u a -V a u p ), r = f |«Y, E p = F ilv u v . 

(3.3) 

These identities follow from the fact that £y = when acting on sources. The vector is 
the acceleration vector, £l^ u is the vorticity tensor, and E^ is the electric field. Since a^, £1^, 
and E^ are all transverse to u^, there are no local gauge-invariant scalars with one derivative. 
Thus, 

Wt = J d d xV=gP(T,v). (3.4) 

We refer the reader to [39] for more details. 

In [19] an alternate method was used to obtain the generating function. By choosing 
coordinates such that V = dt, the d-dimensional metric may be written as 

g = _ e M*)(rfi + a l {x)dx i ) 2 + g^~ 1) (x)dx i dx j , (3.5) 

where gf^ ^ is the metric of the spatial manifold Aid-i- The authors of [19] constructed 
the generating function on Md-i by Wick-rotating to Euclidean signature and then dimen- 
sionally reducing on the time circle. The dimensionally reduced generating function may be 
constructed out of the the dimensionally reduced metric, the gauge field, and the Kaluza- 
Klein (KK) photon aj. For instance, as discussed in section 2.2 of [19], in order to preserve 
gauge invariance, KK invariance, and preserve the form ( |3.5| ) the spatial components of the 
gauge field, A{ must appear in the combination, 4 

Ai = Ai- a t A . (3.6) 

The main difference between the methods of [39] and those of [19] are in the symmetries which 
are made manifest. In [39] the full d-dimensional diffeomorphism and gauge invariances are 
apparent, while in [19] the spatial diffeomorphisms and the symmetry generated by V are 
clearly visible. 

4 In the language of [19] we have A£ urs — A^ heirs , Tours — exp(— , ^ ours = exp(—a)A \ 
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3.2 Chern-Simons contributions to the generating function 

As we have alluded to in the previous section, in addition to scalars built out of gauge-invariant 
tensors, we may add to the generating function terms which are gauge and diffeomorphism- 
invariant up to boundary terms. We presently classify all such terms, as they will play a 
crucial role in the sections to come. Given a conserved, gauge-invariant current X 11 which is 
transverse to we may construct 

W x = J d d xV=9 {ciA^ + c 2 X°) , (3.7) 

where the c% are constants and we choose coordinates so that V = dt- (If X^ is parallel to u M 
then both terms in ( |3.7| ) reduce to scalar expressions.) The first term on the right hand side of 
( |3.7| ) is gauge invariant up to boundary terms. The second term is diffeomorphism- invariant 
up to boundary terms due to the fact that V is Killing. Since Aq is odd under CPT, either 
c\ or c 2 should vanish in a CPT preserving theory. 

In the approach of [19] ( |3.7| ) may be rewritten in the form 

W x = (3 J d d - l x y/g«-V [cxAX 1 ~ c 2 a l X i ) , (3.8) 

where we have defined the Kaluza-Klein invariant covector Xi = e a {Xi — aiXo). Since 
X^u^ = by definition then Xi reduces to Xi = e a Xi. In this notation the spatial indices i,j 



are raised and lowered with the spatial metric ^ defined in (3.5) and its inverse. 



If we denote the Hodge star operator on Md-i by * and treat Ai and Xi as one-forms 
then fl3,8|) can be recast as 

W x = (3 I ciA A *X — c 2 a A *X . 



J c\A A *X — c 2 a A *X . (3.9) 



We would like to argue that *X must either be a product of Chern classes of A and a (dA and 
da) and Pontryagin classes of the Riemann tensor, or that reduces to the gauge-invariant 



scalars and pseudo scalars described in section 3.1 



Conservation of X^ implies that X is conserved, d * X = 0. On a spatial manifold with 
trivial topology there are no non-trivial cohomology groups and every closed form is exact. 
Thus, *X = dy. If y is gauge-invariant then after integrating (|3.9|) by parts we obtain 



W x = P j cxdA Ay- c 2 da A y , (3.10) 

implying that ( |3.8[ ) may be rewritten in terms of the gauge-invariant scalar contributions to 
W described in Section [O]. If y is not gauge invariant, then dy must be gauge invariant in 
which case the gauge variation of y, S\y must be exact on a topologically trivial manifold, 
^a3^ = dZ. Thus, y must a linear combination of d — 3-dimensional Chern-Simons forms 
implying that *X = dy is a product of Chern classes of the U(l) gauge fields and Pontryagin 
classes of the spatial metric. Put differently, the terms in (|3.8| ) are genuine Chern-Simons 
terms on M.d-1- 
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4. Thermodynamics in two dimensions 



In the following sections we will construct the most general Euclidean generating functional 



for a two-dimensional Lorentz- invariant theory to second order in derivatives. In section 144 
we construct the most general generating function to second order in the derivative expansion 
and derive the appropriate constitutive relations. 

We place our theory on a Euclidean cone with metric 

ds 2 = dr 2 + r 2 dr 2 , (4.1) 

where r has periodicity r ~ r + 2tt5 and it parameterizes the Euclidean time circle. The 
parameter 5 is related to the deficit angle of the cone, 2tv(1 — S), thereby interpolating between 
a cylinder (5 — > 0) associated with a thermal state of the theory, and M 2 (5 — > 1) associated 
with the Euclidean vacuum. The expectation value of the stress tensor in the Euclidean 
vacuum vanishes in a conformal theory, and is proportional to the metric in more general 
theories. This property of the stress tensor imposes restrictions on the thermal theory, as we 
now show. 5 

4.1 The generating functional and constitutive relations 

Consider a two dimensional relativistic theory on a manifold Ai with metric g^u in the 
presence of a gauge field A^. We wish to construct the thermal partition function for this 
theory up to second order in the derivative expansion following the algorithm presented in 
Section [| To this end, we need to enumerate all possible scalars, pseudo scalars and conserved 
currents to second order in the derivative expansion assuming that the Lie derivative of the 
metric and gauge field vanish for a timelike vector V*. We have listed all gauge- invariant 
scalars, pseudo scalars and conserved currents to second order in the derivative expansion in 
Table |]. The alert reader may note that no components of the Riemann tensor appear in 
Table |l[ This is because of the fact that, in a two-dimensional time-independent background, 
the Riemann tensor is given by 

R^vpa = {J; ~ si) (svsw - g^9up) , (4.2) 

where s\ and S5 are defined in Table [l], and so the Riemann tensor is not independent. Thus, 
to second order in derivatives the most general generating functional takes the form 



J d 2 x^~g (p(T, pL) + c 2d (3- 1 X° + sidi(T, //) + ^ s iai (T, /z) j . (4.3) 



Wo 



We have not included the scalars S4 and S5 in ( [4. 3D since they are total derivatives which, after 
partial integration, contribute to ol\ and 03. Similarly, there exists a linear combination of 



5 The S — > limit seems somewhat singular in the coordinate system ( |4.l[ ). This can be remedied by going 
to a coordinate system defined such that r — j-e^f and 8 = <j>8 which is valid for any non zero 8 where the 
line element takes the form ds 2 — e~fT (dp 1 + R 2 d<f> 2 ) whose 5 — > limit gives us the cylinder. 
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si and §2 which vanishes after integrating by parts. Thus, we have omitted a term involving 
%. The coefficient of the Chern-Simons term X° has been denoted C2d to distinguish it from 
the coefficient of the Chern-Simons term of a four dimensional theory to be introduced in the 
next section. 



1st order 


2nd order 


conserved currents 


h = 


si = a^a p 
s 2 = F^F^ 
s 3 = V„TW# 
S4 = V„V"£ 
s 5 = V^WT 


X p = Tu p 



Table 1: A list of all gauge- invariant scalars and conserved currents formed out of the background 
fields and V with at most two derivatives. The quantities u^, T, \l and a M were defined in Section 3.1. 
The transverse vector is defined through £t M = e^ v u v . 



If the theory has a U(l) or a gravitational anomaly then we must add to the generating 
function extra terms which break gauge and diffeomorphism invariance such that the resulting 
current and stress tensor are anomalous. We find that the appropriate terms we need to add 
to the generating function are given by 

W A = J d 2 aV=S (-c s Mm«" - c 9 u Q ?A^rv) • ( 4 - 4 ) 

In [23] the second expression in ( [4.4 ) has been obtained by integrating the equations of 
motion. The covariant energy momentum tensor Tc u v and current J cov are given by ( [2.15 ) 
with W = W2 + Wa, and one can check that the resulting (non) conservation equations are 
given by Q2.37D , 

V — r f^ v F V T pv — F v T p A- r e up 8 /? (A 51 

v H u cov ~ t -« t r pi v fJl ± cov — r p'J C ov' c 9 t u p ix. K^- ) 

An explicit computation gives the one-point functions of the covariant current and stress 
tensor to be 

J p ov = R u p + Ru p , TZ = Eu p u u + Vu p u u + 9 {u p u u + u v u") (4.6) 
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where we define was denned in the caption of table |] and 

dP dai_ (dai da 3 da 3 \ da 2 ( da 2 da 2 \ a 3 

R = - 2c s ^, 

e v , dP , rr^dP ldai, ( da, ^daA ( da 2 da 2 1 da 3 \ 

da,\ 2a, 
as + 2-^- I s 3 - Ta 3 s 4 — s 5 , (4.7) 

=P — a,s, — a 2 s 2 — a 3 s 3 , 
6 =~c 2d T 2 - c sf i 2 + ^s 5 . 



4.2 Conformal theories 

In a conformal theory and in the absence of anomalies, the generating function must be invari- 
ant under Weyl rescalings of the metric. In terms of the parameters in (^^) we find that this 
implies that P = T 2 po(n/T), a\ = 0, 6 a± = 03 = and that a<i = a2{n/T)/T 2 . In addition 
to Weyl-invariance, conformal symmetry enlarges the U(l) gauge symmetry to a doublet of 
current algebras for a left-moving and right-moving current. This doublet gives additional 
constraints on the generating functional, owing to the independent conservation laws for the 
chiral currents. For instance, the chiral anomaly fixes the dynamics of U(l) currents when in 
the high temperature limit, when the chemical potential vanishes, or when the £7(1) is non- 
compact [31,51]. (Compact U(l)s may be more interesting due to non-perturbative effects, 
like the twisting of winding modes in the presence of a chemical potential.) In the remainder 
of this section we will ignore the effects of the doublet of current algebras and leave them for 
future inquiry. 

In the absence of a U(l) current we would have a<i = so that all first and second order 
terms in the derivative expansion are excluded. In fact, conformal invariance precludes us 
from writing any higher derivative correction to the generating function. Indeed, suppose that 
there exists a scalar expression which contributes to the generating function. Such a scalar 
must be constructed out of the Weyl tensor W^ upa , the velocity field u^, the temperature T 
and their gradients. Since gradients of Weyl covariant quantities generate non Weyl covariant 
expressions it is convenient to trade the gradient with a Weyl covariant one Vj^ defined in 
[52]. Since V^u u = 0, Vj^T = and the Weyl tensor identically vanishes in two dimensions, 
there are no Weyl-covariant tensors with one or more derivatives. A moment's thought also 
shows that we may not add any terms non-analytic in derivatives either. Furthermore, the 
analysis of section |3.2| implies that there are no conserved currents other than X^ 1 defined in 

6 In practice, if we set cti = &i(T) then the expression — gai§i becomes a total derivative and does not 
contribute to the energy momentum tensor and current. Similarly, had we used 02(/x/T)s2 in place of ai(T)§i 
we would have found that it contributes only to a boundary term. 



- 21 - 



table |I|. It then follows that the equilibrium partition function of an uncharged conformal 
fluid contains only the pressure and anomaly terms. 

Often, the Weyl symmetry of a conformal theory will be anomalous. This anomaly 
manifests itself as a non- vanishing trace of the stress tensor 

(T co ,X = c w R, (4.8) 



where c w is proportional to the central charge as defined in (1.9) and R is the Ricci scalar. 
The Weyl anomaly Q4.&] ) can be generated by allowing for a non- Weyl covariant contribution 
to the generating function of the form a\ = c w . Thus, in the presence of gravitational, gauge 
and Weyl anomalies, the generating functional of a conformal theory will take the form 

W 2 = J d 2 x^~g (r 2 p (/i/T) + hd^X + a2 ^ )s2 + CwSl ^j + W A , (4.9) 



where Wa was given in ( [4.4|) . Using (4.7) we can write the resulting energy momentum tensor 
and current in the form 

p dP a' 2 4a 2 
R = - 2c s/ u, 

£ =poT 2 + c w si - ||s 2 - ^«5, (4.10) 
V =pqT - c w si - j^s 2 , 
=c 2d T 2 - c s fi 2 + ^-s 5 . 

Note that 

V M V^ = -^d^sTdu^ (4-11) 

transforms homogeneously under Weyl transformations with weight +2 in two dimensions. 
As explained previously, in the absence of a U(l) current the relation ( [4.10 ) is exact. We note 



that relation (|4.10| ) is exact also in the absence of an external electric field wherein gradients 
of the chemical potential are proportional to the acceleration such that V a (fJ./T) = (and 
then S2 = ^4 = 0). 

Let us place our theory on a cone with line element 



ds z = dr 2 + r z dr 2 , (4.12) 

with r ~ t + 2tt5 so that 27r(l — 5) is the deficit angle of the cone. The limit 5 — > 1 corresponds 
to Euclidean space and (5—^0 corresponds to the cylinder. Treating the angular coordinate 
as the Euclidean time direction, the temperature of the theory on the cone is given by 

T" 1 = 2ir5r . (4.13) 
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Thus, we have, for example, V^T _1 V At T _1 = 4-7r 2 <5 2 , and the gradient expansion on the cone 
breaks down. Despite the breakdown of the derivative expansion, the theory on the cone is 
in a type of time-independent equilibrium state. In this time-independent state, we continue 
to formally define the temperature T as the inverse length of the time circle. As we have just 
argued, the generating function W is the most general Euclidean generating function on can 
write down for a two dimensional conformal theory. Therefore, its variation will correctly 
capture the expectation value and current at least far from the tip of the cone. 
For an uncharged conformal theory on the cone, 

TZ = iPo ~ 4n 2 c w 5 2 )T\2u^ + g^) + (c 2d + 8ir 2 5 2 c g )T 2 (u^ + u v u») . (4.14) 

In the limit where 5 = 1 we expect to recover the energy momentum tensor in the Euclidean 
vacuum of a conformal theory on R 2 which should vanish. Thus, 

po = 4:TT 2 c w , (4.15a) 
c 2d = -8ir 2 c g . (4.15b) 

We note in passing that instead of placing the theory on a cone, one could have placed the 
theory on a punctured plane M 2 /{0, 0} and carried out the entire calculation by setting 5 = 1 
from the outset. One would then require that after adding the missing point at the origin, 
the stress tensor ( |4.15 ) should coincide with the stress tensor of the Euclidean vacuum, i.e., it 



should vanish. We make some further comments about this alternate viewpoint in section 6.2 . 

Once there is a £7(1) current, we need to account for higher derivative terms in the 
generating functional and one might worry that these will spoil ( [4.15b| ). However, we are 
interested only in those terms which do not vanish when we set the external gauge field to zero 
so the only worrisome terms which we need to account for are those that contain gradients 
of the chemical potential. As noted previously, in the absence of an external electric field 



V a (///T) = 0. Thus, ( 4.15a ) remains valid to leading order in the chemical potential, and 



( [4.1 5b| ) does not get corrected 

po = Att 2 c w + 0{n) , c 2 d = -8ir 2 c g . 
4.3 Non-conformal theories 

An important ingredient in obtaining ( |4.15| ) was the absence of higher derivative contributions 
to the generating function. Had there been higher derivative corrections to the generating 
function, the breakdown of the derivative expansion on the cone would have forced us to 
consider the possible contribution of any such corrections to the energy momentum tensor. 
We now turn to such considerations. 



Using (4.7) to compute T^ v explicitly, we find 



c 2d + 8ir 2 5 2 c g 
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We would like to argue that T^ v does not receive corrections at any finite order in the 
derivative expansion. To this end, consider the perturbed cone, given by 

g = dr 2 + r 2 dr 2 + h Tr (r)drd,T . (4.17) 

That corrections to T^ v are absent amounts to the statement that there are no scalars or 
conserved currents which we can add to the generating function which are linear in h rr . We 
will now prove such a statement. 

In the absence of external gauge fields, the tensors that may contribute to the Euclidean 
generating functional are built out of u M , T, /i, R^upa, and derivatives thereof. In two dimen- 
sions the vorticity tensor vanishes so that the general expression (|3.2[) becomes 



V^T = -Tap, V M — = V^Uj, = -u^a v . (4.18) 



The Riemann tensor in two dimensions takes the form Q4.2p , and on the cone, to linear order 
in h Tr , 

V^a^ = -u^u v a 2 - a^a v . (4-19) 

Thus, the most general scalar expression in a two dimensional thermodynamic theory may 
be constructed out of a^, u^, the metric and the epsilon tensor e^ u (and no derivatives). 
A quick computation then shows that one can not construct a scalar linear in h Tr . Following 



the discussion in section 3.2 there are no conserved currents in two dimensions other than 
which will contribute non trivially to the generating function. Therefore, at any finite order 
in the derivative expansion there are no gauge or diffeomorphism invariant expressions which 
can contribute to T TV on the cone in the absence of sources. We also note that our argument 
precludes local non perturbative terms from contributing to T rr as well. For example, an 
expression of the form e~~ , with T the temperature s a scalar and n some number, can not 
contribute to T TV since s has no contribution linear in h Tr when evaluated on the background 

(£!)• 

At 5 = 1, T^ v should vanish by the translational invariance of the Euclidean vacuum. 
Since T^ v is given precisely by ( [4.16 ) then requiring it to vanish in the limit 6 —> 1 implies 
that 

C2d = -8-K 2 c g . 

This demonstrates our claim that ( [4.15b| ) remains valid for non conformal theories. 
5. Four dimensional theories 

We move our attention to four dimensional theories with anomalies. We start by constructing 
the gauge and diffeomorphism invariant generating function to first order in derivatives. As 



mentioned in Section 3.1, there are no local gauge- invariant scalars with one derivative but 



there is one (and only one) CPT preserving Chern-Simons term of the type described in 



section 3.2 given by a A dA. The corresponding conserved current is 

X? = T(B» + ^), (5.1) 
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where the magnetic field B p and vorticity are defined by 



B" = -e^u u F pa , = e^u v V p u a , 



(5.2) 



In the absence of anomalies, the generating function to first order in derivatives is given by 

Wi = J d 4 x^9 {P(T, n) + c u r X Xl) . (5.3) 

The alert reader will note that there exists an additional conserved current X p = T 2 uj p which 
arguably contributes to the generating function a term of the form Wx = f d^x^—gcA^X^ ■ 
However, the scalar A p X% corresponds to the Chern-Simons term A A da which is equivalent 
to a A dA after integrating by parts. Indeed, f d A x^—g(3~ 1 X® = f trxy '—gA^Xrf up to 
boundary terms. 

If the theory has t/ (l) 3 and mixed anomalies, the partition function has an anomalous 
variation given by (minus) ( p. 41 ). Such an anomalous variation is generated by 



W A = j d A x^f^A p (c A f A + Cmijfa - ajcs)) 



(5.4) 



with 



Ja 



Jm 



-2// {B^ + ^ 



Pes = ^ 



i^W pvpa u v u p uj a + X - (r + 6R up u u u p - 6a 2 - ^lu 2 ^J oj p 
v v ft r T 4- — r v r T I 

1 VT U P L (TV ^ VT p<f> L (TV J ' 



(5.5) 



where W pvpa is the Weyl tensor. The covariant current and stress tensor can be obtained by 
varying the generating function W = W\ + Wa and adding the appropriate Bardeen-Zumino 
polynomials ( |2.45 ), as defined in ( [2.15 ). The stress tensor and current obtained this way 
satisfy the anomalous conservation equations (2.46), 



1 



M con 



3c A F Kcr F a {3 -\- C m R \kctR va/3 



e™ ap F Ka R» v a p 



(5.6) 



Y7 7^ — zpii ju _i_2r V 
The resulting constitutive relations for the stress tensor and current are presented in appendix 



Consider the TIL, component of the stress tensor on a cone times M 2 , 



ds 2 = dr 2 + r 2 dr 2 + dx 2 + dy 2 . 



(5.7) 



where r ~ r + 2tt5, with a magnetic flux given by the gauge potential A y = Bx with B 
constant. We would like to argue that, in this background, the T rr component of the stress 
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tensor can only receive contributions from the terms present in (5.3) and fl5.4|). To see this 



we turn on a metric perturbation 5g Tr = h Tr (r,x,y) and look for terms in the generating 



function which are linear in the perturbation. Recall that, following section 3.2, we have 
already accounted for all possible Chern-Simons contributions to the generating function W 
allowing us to focus our attention on gauge-invariant contributions only. In practice, since 
all gauge-invariant variables do not depend on the M. 2 coordinates x and y, it is sufficient to 
consider 5g Tr = h Tr (r). Indeed, all spatial derivatives of the M 2 coordinates may act only 
on h rr and will therefore lead to a vanishing contribution to the generating function after 
integrating by parts. For instance, consider the three derivative contribution to the generating 
function of the form 



W 



J (fixyf^R^uVa? = ip J drdxdy (d 2 x + d 2 ) ^ . (5.8) 



Thus, 6W/5g Tr (r,x,y) = 5W/5g Tr (r). 

Focusing on perturbations of the form 5g Tr = h rr (r) the Riemann tensor, the vorticity 
tensor Q^u an d the electric field vanish at least to 0(h 2 ) and 

V M u„ = -u^dy + 0(h 2 ) , 
V U T = —Ta a , 

2 ,2, ( 5 - 9 ) 

V = -u^UyCL - a^a v + 0{h ) , 

V^By = - a^B u - a a B a u fl u u + 0(h 2 ) . 

It is a straightforward procedure to check that there are no scalars or pseudoscalars which 
can be constructed out of u^, B^ and a M which are linear in h Tr . As was the case in two 
dimensions, the absence of scalars and pseudo scalars which are linear in h Tr implies that 
local non perturbative contributions to W will not affect TJ| 



^COV 



We have argued that in the background given by (|5.7|), the only contribution to T c 



- cov 



comes from ( |5.3D and (|5.4|) . An explicit computation (see |B.3|) gives us 



T cov = * B 4n 2 § 2 r 3 ' ( 5 - 10 ) 

We then use the same argument as in the two dimensional case, namely that the translational 
invariance of the Euclidean vacuum implies that T^ v vanishes in the limit S — > 1, to obtain 

Cid = -8ir 2 c m . (5.11) 



6. Summary and discussion 

Thermodynamic response coefficients induced by gravitational anomalies pose entirely new 
types of challenges when compared to response coefficients induced by pure U (1) anomalies. In 
particular, the gravitational anomaly coefficient c g in two dimensions and the mixed anomaly 
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coefficient c m in four dimensions were conjectured to be related to the respective response 
parameters C2d and via 

c 2d + 8vr 2 c 9 = (6.1a) 
c M + 8vr 2 c m = . (6.1b) 

These conjectures have a natural generalisation to any even dimension [29]. They are peculiar 
in that (a) they mix coefficients of different orders in the derivative expansion and (b) they 
involve the transcendental number n. A main aim of this work is to fill up this lacuna in our 
understanding of gravitational anomaly induced response and to prove (|6.1|). We will begin 
our discussion in section [O] by reviewing our derivation of ( |6.1] ). Then, in |6.2j we critically 
examine various assumptions associated with our technique. We end with a discussion of 



future prospects in section 3.2 



6.1 Thermodynamics on cones and gravitational anomalies 

Consider a two dimensional Euclidean theory on a cone C$ of deficit angle 2tt(5 — 1). The 
Euclidean partition function, W$, of such a theory can be thought of as the thermal partition 
function (at a temperature T = (27r<5r) _1 ) with the angular direction being the Euclidean 
time. Since \d r T\ = (2it5)T 2 rather than the required parametric suppression for the deriva- 
tive expansion \d r T\ -C T 2 , the derivative expansion breaks down on the cone. (This break- 
down of the derivative expansion is unrelated to various problems that may originate from 
the conical singularity at the tip of the cone, which we will return to in the next subsection.) 

A priori, the failure of the derivative expansion means that one cannot directly use 
the methods of [19, 39] which parametrize thermal partition functions on arbitrary time- 
independent backgrounds using the derivative expansion. Despite this, we argue that one 
may use the methods in [19, 39] to accurately compute a specific off-diagonal component in 
the one-point function of the stress tensor T TV of the Euclidean theory. Further, we argue 
that T rr receives contributions from only two terms among the infinite number of terms in 
the derivative expansion including any local non perturbative contributions to W. The two 
terms on which T rr depend are the contribution of the gravitational anomaly proportional 
to c g and a certain Chern-Simons term associated with C2d- 

When there is no deficit angle our generating function reduces to that of the Euclidean 
vacuum, Wy ac , which is presumably unique. Assuming that the generating function is con- 
tinuous at 5 = 1, 

W Vac = lim W 5 , (6.2) 
<5— >1 



we demand that lim^i T rr = resulting in ( |6.1a ). One way of interpreting our computation 



is to think of T rr as a Casimir momentum of the thermal state which is generated in taking the 
Euclidean theory from the plane R 2 to the cylinder S 1 x R. Phrased this way, our argument 
is conceptually similar to the argument used to derive the Cardy formula. Indeed, in the 



framework of conformal field theories we can prove ( 6.1a ) using a formalism similar to that 
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of [40, 41] — see appendix [A] for details. We emphasize that, unlike the Cardy formula, our 
argument involving the generating function can be applied to non-conformal theories. 

Our derivation of ( |6. la|) (and the Cardy formula argument) satisfactorily circumvent 
both of the a priori obstacles mentioned in the introduction faced by previous methods: the 
transcendental factor of 8ir 2 enters into the argument geometrically through the radius of 
the Euclidean circle, and the breakdown of the derivative expansion on the cone allows the 
gravitational anomaly coefficient c g which is second order in the derivative expansion and the 
coefficient of the Chern-Simons term, C2d which is zero order in the derivative expansion to 
contribute to correlation functions at the same order in derivatives generating a "jump" in 
the derivative expansion. 

The situation in four dimensions is very similar to the two dimensional case. To derive 
the constraint (6.1b) involving the mixed gravitational anomaly coefficient we place our the 



theory on cone times a plane, x M 2 , with a small but a constant magnetic flux turned on 
over M 2 . We then consider the off-diagonal component of the stress tensor whose two indices 
are on the cone. As in the two dimensional case, we can argue that only the terms associated 
with did and c m terms contribute to T rr . Using the same assumptions as those spelt out for 
the two dimensional case we arrive at ( |6.1b ). 



6.2 Validity of our results 

In this subsection we would like to examine possible pitfalls and the validity of various as- 
sumptions underlying our derivation of fl6.l| ) . The starting point of our analysis was to place 
the field theory under consideration in a background with a conical singularity. To do this in 
practice one must usually regulate the singularity by smoothing out the conical tip in some 
particular way. Following [53] one may consider a class of regulators of the form 

dsf = r 2 dr 2 + F{r/l)dr 2 (6.3) 

where Fix) is a smooth function of its argument and satisfies F(0) = 5 2 and F(oo) = 1. 
For instance, F = (r 2 /I 2 + 5 2 )/(r 2 /£ 2 + 1) which describes a hyperbolic space satisfies these 
criteria. 

Since the sequence of geometries ( |6.3| ) are smooth they can not converge uniformly to the 
cone which is singular. Indeed, one of the main results of [53] is that while various geometric 
observables (such as arbitrary polynomials of the curvature tensor) depend on the details 
of the regularization procedure, there are other observables which are independent of such 
details. We expect this to be true of field theory observables as well. In particular, we expect 
to find some arbitrariness in the structure of possible localized states near the singularity. 
(This phenomenon is well-known from string theory computations involving resolution of 
orbifold singularities.) 

Whether they depend on the regularization parameter or not, the distressing feature of 
localized states at the tip of the cone is that the Euclidean generating function on the cone, 
Wg, will not approach the vacuum generating function in the limit 5 — > 1, i.e., Wg is not 
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continuous at 5 = 1. Since our analysis relied on taking the 8 — )■ 1 limit of the generating 
function as in (|6.2|), one might worry that our computation would fail in the presence of 
localized states at the tip. However, our analysis focused on computing the T TT components 
of the stress tensor far from the tip of the cone. Therefore, if the equality ( |6.2| ) is true up 
to localized expressions at the tip then our computation remains unspoiled. In other words, 
since we are computing components of the stress tensor away from the tip of the cone, we may 
require a weaker condition than (S.2). Namely, that the equality in ( |6.2| ) holds up to terms 
localized at the tip. Thus, our computation is robust against the addition of any localized 
state near the singularity. 

Despite our computation being robust under the addition of any localized state near the 
tip, one should still be concerned about the continuity of W$ at 5 = 1. We emphasise that 
this is not an empty concern. As is known from heat kernel computations on the cone [54] the 
heat kernel (and hence the partition function) of massless higher spin fields (gravitinos and 
gravitons in this case) are not continuous at 5 = 1. The reason for this discontinuity is the 
enhanced symmetry on the plane. On the cone we have only one Killing vector, the generator 
of rotations, while on the plane we have three Killing vectors which generate ISO(l, 1), two 
translations and rotations. Since the graviton and gravitino are susceptible to the geometry 
the number of zero modes of these fields is discontinuous at S = 1. This suggests, for 
example, that the relations (^0|) will fail for the gravitational anomaly contribution from 
chiral gravitino. This was in fact known before from explicit computations of C2 n d f° r a f ree 
theory of chiral gravitinos [55]. It is satisfying that our derivation excludes those systems 
where the relations in question are known to break down. 

Instead of placing the theory on a cone Cs we could have placed our theory on M 2 '* = 
M 2 /{0, 0}, i.e., the plane with the origin removed. In this case, instead of fl6.2| ) we would have 
had to assume that the origin can be added back without modifying the generating function 
W away from {0,0}. In other words, If the generating function W on the punctured plane 
differs from the generating function W on M? only by terms which are localized at the origin 
then our analysis carries through without any obstruction. 

This discussion leads us to conclude that the relations (|6.l|) hold for all local field theories 
lacking massless gravitini or gravitons. This is certainly true for free field theories [54] and 
therefore also for weakly coupled ones. It would be satisfying to relate our argument to a 
mathematical theorem for a specific index which would bypass this problem of continuity 
entirely. 

As we have emphasized several times, the derivative expansion breaks down on the cone. 
As a result, it is insufficient to use a truncated generating function W m valid only up to order 
m in the derivative expansion. Let us restate how this problem has been circumvented. Owing 
to the fact that in Euclidean space correlators drop of exponentially at small distances as in 
(3.1), the full Euclidean generating function may be formally expanded in a power series in 
derivatives around the flat space thermodynamic equilibrium configuration. When resummed, 
this infinite series together with any non analytic terms at the origin will reconstruct the 
Euclidean partition function. What we have argued for in sections |] and [5] is that there are 
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no terms in this series and no local non analytic terms at the origin which may contribute to 
the particular component of the stress tensor which we were interested in. 

6.3 Future prospects 

We will now conclude with various thoughts on future work. We begin by noting that the 
parameters c^d and c^d that are the focus of this paper are respectively the coefficients of the 
Id pure and 3d mixed Chern-Simons terms for the thermal Kaluza-Klein field. The relations 
( |6.lD essentially relate these Chern-Simons terms to pure or mixed gravitational anomaly 
coefficient in the theory with one dimension higher. 

A relation between Chern-Simons terms and anomalies is already familiar in the case of 
U(l) anomalies [56] where an abelian anomaly in 2n dimensions is related to a parity anomaly 
in 2n — 1 dimensions. It would be interesting to see whether a similar argument could be 
used to directly derive Kaluza-Klein Chern-Simons terms in field theory after compactifying 
on a warped thermal circle. We refer the reader to [42] for related discussions. 

Next, the procedure that we have used to establish ( |6.1| ) relied on properties of the 
partition function on a cone in the limit where the deficit angle goes to zero. This object has 
been studied in the literature on entanglement entropy (see [57] for a review). Briefly, given 
a state in, say, a 3 + 1 dimensional field theory, and a plane which cuts the three dimensional 
space into two halves, the entanglement entropy of half the space is the Von-Neumann entropy 
associated with the density matrix obtained after integrating over the of degrees of freedom on 
one half of the space. An efficient way to compute this entanglement entropy is via the replica 
trick (which for a field theory leads to the so called "conical singularity" method) whereby 
one introduces a conical deficit on the interface and studies the resulting partition function 
as a function of the deficit angle. It would be interesting to see whether our analysis can be 
translated into specific statements about how gravitational anomaly enters the entanglement 
entropy. 

While we have argued for ([O]) which are valid in two and four dimensions it would be 
interesting to try and prove the higher dimensional generalisation of ( |6.lD proposed in [29]. 
There, a "replacement rule" has been proposed where, in the current language, the leading 
order Chern-Simons coefficients in the generating function are related to Chern-Simons terms 
and first Pontryagin classes of the anomaly polynomial. It seems more than likely that one 
may be able to prove the replacement rule by placing a 2n dimensional theory on spaces with 
particular fluxes effectively reducing the problem to the two dimensional one which we have 
solved. Furthermore, one can hope to generalize the replacement rule of [29] in such a way 
that higher order Pontryagin classes in the anomaly polynomial relate to Chern-Simons term 
in the generating function which are higher order in derivatives. We leave these and other 
related problems for future work. 
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A. Casimir energy density and energy flux in two dimensional conformal 
field theories 



Consider a conformal theory with central charges cl and cr on an arbitrary manifold with 
line element 

ds 2 = p(z, z)dzdz . (A.l) 
Using R = — Dlnp, the solution to ( |2.37| ) is given by 

rpZZ „ „— 1 



c w p Dlnp, 
TTL = 2{Cw ± 2Cg) [2pdlp - 3(d z p) 2 ) , 



rpZZ 

1 cov 



p 



{2pd 2 zP -3{d z p) 2 ) , 



(A.2) 



where we have required that Tcov = for p = 1 and defined the linear combinations 



Con -\~ *2Cq 



CL 

24vr 



2c„ 



cr 

24vr ' 



(A.3) 



in accordance with the conventions of the main text, namely equations ( |2.28| ) and (|1.9|). 

A cylinder of radius R can be parameterized by the line element ds 2 = dx 2 + d9 2 with 
9 ~ 9 + 2tt/R which can be obtained from (A.l) using 



p=—, z = Re { - x+i9 ^ R . 

zz 

Inserting (A.4) into ( |A.2j ) and Wick rotating to Lorentzian signature we obtain 

47T 2 c„, — 8ir 2 c„ \ „ 9 



(A.4) 



COV 



-M 2 Cg 4:TT 2 C n 



(A.5) 
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with T^ 1 = 2vri?. Thus, the energy momentum tensor of a conformal theory on a cylinder dif- 
fers from that of flat space by a Casimir energy and Casimir momentum density proportional 
to combinations of the left and right moving central charges. Therefore, the energy density 
and energy flux of a conformal theory are completely specified by the central charges of the 
theory even though these enter at second order in the derivative expansion. The result flSjgp 



generalizes the results of [40, 41] to theories with a gravitational anomaly; the off diagonal 
term of the energy momentum tensor has a universal value on the infinite strip proportional 
to the difference between the left and right moving central charges. 

B. Constitutive relations associated with anomalies in four dimensional the- 
ories. 

In this section we list the part of the four dimensional constitutive relations which follow from 
the anomalous part of the generating functional, Wa, which includes the three derivative terms 
proportional to the mixed anomaly coefficient c m . This is the four dimensional analogue of 
the two derivative terms derived in [23] for two dimensional theories with a pure gravitational 
anomaly. 

We will now focus on contributions to the covariant stress tensor and current obtained 
from Wa in equation (5.4) by variation, according to the definition ( |2.15| ). We will denote 



them by T^ cov ^ and Ja CO v respectively. We will present these contributions in the following 
form 

with qAi^A-, t~a transverse and ta traceleless, and the subscript A denotes the contribution due 
to the anomalies. We find that the resulting expressions for £a, Va, v% Q% and contain 
a number of Weyl-covariant terms, as well as a number of Weyl-non-covariant terms. 7 For 
this reason, we have found it useful to summarize a number of independent Weyl-covariant 
tensor structures in Tables |2| and |3|. There are also the Weyl-covariant zeroth order quantities 
T, n and , as well as the first order pseudovectors and defined through 

uP = e^UvVpUv, B* = X -e^u v F pa , (B.2) 

and the antisymmetric tensor £l^ v defined in the caption of Table |2[ 

In terms of these Weyl covariant quantities, we find that the scalars Ma, £a, and Va are 
given by 

Ma = 0, £ A = 3V a = 2c m w„(f«# 1 - w£ 2 ) + 4c m B li {v^ 1 - W v a v ). (B.3a) 



7 Because the gauge and mixed gravitational anomalies may be non- vanishing in a CFT, we find it interesting 



that there are non- Weyl-covariant terms in (B.l). We leave a comprehensive understanding of these terms and 



their relation to the Weyl anomaly for future work. 
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1 


2 


2nd order scalar (s2,i ) 
2nd order vector {v 2 %) 
2nd order tensor (t 2 ) 


R + 6R pu u p u u - 6a 2 
A^R up uP - 2VL pv a v 
W pp(jv u p u a 


A^V v F pv 



Table 2: Some independent Weyl-covariant second-order tensors, where we have defined antisymmet- 
ric tensor Q pL/ = A /ip A 1/0 '(VpU (T — W a u p )/2. The other independent second-order tensors are products 
of first-order tensors. 







3rd order pseudovector (S3) 
3rd order pseudotensor (£3) 


e pupr7 u u (V p Raa)u a + a^ v AP^V u u p + u p a 2 
A *<n e v>P*T Uf) iy aRra + 2a a t 2Ta ) 



Table 3: Some Weyl-covariant third-order pseudotensors. 



The vectors v A and q A are given by 

- 6c Af xBi" - 3caH 2 ^ - Ac m t^u u - ^c m (s 2A - ^ooA uj p , (B.3b) 



ft 
A 



q\ = - 3c Af i 2 B p - 2c A fi 3 u p - 2c m U^B V + - a 2 ) B^j 

- 2/ic m (vg + ttfuv + i (s 2A - 00A + c m {A pv p u v + 2u p E u a v - E p a v u v ) , 

where we have defined the transverse tensor § pv = R ppua u p u a . Finally the tensor is 

t% =4^c m (if - W^^u^ + 2oj <p v^>) + 2c m w < ' i u0 + 4c m uj <p n u>p E p (B.3c) 
+ 4c m B<^0 + 2c m e^ < %£ (7 (2^ > - fl^A"* 1 ) 
- 2c m V < »F v>p u p + 4c m a < ^ > ^ pJ B (T u; T , 

where the angular brackets denote the transverse traceless projection 

v <fiv> = A w A ""^ ff) _ -A p,/ A pa V pa . (B.4) 

3 

As expected from the arguments given in section |l[ all of the expressions in ( |B.3| ) proportional 
to c m have three derivatives and so do not relate the mixed anomaly coefficient c m to the first 
derivative term listed in J1.6p . 
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